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I. INTRODUCTION 


A. DESCRIPTION OF THE TRANSOGNIC COMPRESSOR TEST RIG 

The Transonic Compressor Test Rig at the Turbopropulsion 
Laboratory (TPL) of the Naval Postgraduate School is shown 
schematically in Fig. 1 and consists of the following major 
components: 

1. Air drive turbine. 

2. Air supply system. 

3. Associated piping including throttling valves at the 
turbine and compressor inlets. 

4. Test compressor. 

The drive turbine is a dual-flow axial air turbine with 
50% reaction. The geometry is given in Table 1. The profile 
shapes of the turbine rotor and of the stator blades are 
identical and the blades are of constant section along the 
radius as shown in Fig. 2. The stator has 31 blades while 
the rotor has 32 (to avoid resonant excitation from wake inter- 
ference). The two parallel stages of the turbine are designed 
for the following output and total inlet conditions: 

Pressure Ratio: 2.8 

Total Inlet Temperature: 640°R 

Flow rate: 10.85 LBM/SEC 


Horsepower: 485 HP 
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The compressor presently under test is a transonic single 
stage, axial flow compressor. It is instrumented for measure- 
ments of torque, mass flow rate, stagnation temperatures and 
pressures, case and hub wall pressures, and for unsteady 
pressure measurements in the flow field and at the walls. 

The Air Supply System incorporates an electric motor- 
driven multi-stage axial flow compressor manufactured by 
Allis-Chalmers. It can presently supply up to 12 lbs/sec of 
air at 3 atmospheres, at temperatures between 560°R and 660°R. 
The compressor is rated at 1250 HP and has a controlled 


variable speed drive. 


B. STATEMENT OF THE TASK 

The Transonic Compressor Test Rig was designed to pro- 
vide the means for obtaining experimental data in fundamental 
compressor phenomena. Following the present experiments, an 
experiment to investigate the onset of supersonic unstalled 
blade flutter is planned which would involve replacing at 
least the present compressor rotor by a rotating cascade of 
flat-plate blades. Such a rotor would not be able to produce 
the pressure ratios required to pump the required flow rates 
through the system. Therefore, it has been proposed, that 
a turbocharger compressor be fitted in series with the 
rotating cascade to provide the required flow through it. 
The turbocharger would also be driven using air from the 


Allis-Charmers air supply system. 
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In order to evaluate the feasability of the turbocharger 
installation, it is necessary to determine the mass flow rate 
required by the drive turbine to drive the test compressor 
at a given power and speed. The remaining air to drive the 
turbocharger turbine is then known and the selection of a 
commercially available turbocharger suitable for this 
application can be made. 

Thus, the performance of the air drive turbine must be 
known over the complete speed range. Of particular impor- 
tance, are the required mass flow rates for given values of 
horsepower. The problem, therefore, is to obtain the turbine 


performance map for all pressure ratios and speeds. 
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II. APPROACH 


A. BACKGROUND 

A search of the most recent literature revealed a number 
of analytical methods for the calculation of turbine off- 
design performance. The majority of these used in a finite 
element approach but little information on the relative suc- 
cess of these methods in practice was available. Two alternate 
methods, both used at the Turbopropulsion Laboratory, were 
those of M. H. Vavra and E. Macchi. Each was examined in 
detail. 

The method of Vavra, given in Ref. [1] is a one-dimen- 
Sional (meanline) approach using mathematical modelling and 
experimental data to express flow angles and losses. It is 
primarily a method to design turbine blading but may also be 
used to predict turbine performance for a given set of gas 
inlet and operating conditions when the blading geometries 
are specified. It is assumed that the axial velocity is 
constant along the blading from hub to tip. Vavra states 
that this assumtion is reasonable for blading in which the 
tip-to-hub ratio is equal to or less than 1.15. The ratio 
is 1.312 and 1.424 for the drive turbine stator and rotor 
blading respectively. It was thought therefore, that the 


method of Macchi might yield more accurate predictions. 
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Macchi's method is given in Ref. [2]. The method, imple- 
mented by Macchi in a computer program written for the IBM 360, 
was an extension of the work done by R. Eckert [Ref. 3] and 
R. Harrison [Ref. 4}. Eckert wrote a program, following a 
simplified three-dimensional analysis, which could be used to 
predict the performance of a single-stage axial flow turbine. 
Harrison improved the program by modifying the analysis to 
take into account streamline curvature. Both programs were 
based on the three-dimensional method developed by Vavra in 
Ref. [5]. Macchi'ts principle improvements to the program were 
to introduce the choice of various methods to calculate gas 
outlet angles and loss coefficients. Two methods of calcula- 
ting gas outlet angles are included; those of Ainley and 
Mathieson [Ref. 6] and Traupel [Ref. 7]. Five methods for 
calculating the loss coefficients can be selected; those due 
to Ainley and Mathieson [Ref. 6], Dunham and Came [Ref. 8], 
Balje [Ref. 9], Lonherr and Carter [Ref. 10] and Traupel 
{Ref. 7]. 

Macchi's computer program, as documented in Ref. ]2], was 
selected for performance predictions of the drive turbine. 

It should be noted that no card deck of the program was avail- 
able, and no results of using the program were available 


other than those included in Ref. [2]. 


B. ANALYSIS 

The method requires the following assumptions; 

1. There are an infinite number of blades in each blade 
row so that blades downstream do not affect upstream conditions. 
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2. The flow is axisymmetric at locations where the 
equation of motion is solved. 

3. The flow is steady and adiabatic. Thus, the total 
enthalpy through the stator remains constant along a stream- 
line and the relative total enthalpy through the rotor remains 
constant along a streamline. 

4. All equations are solved at between blade row loca- 
tions. Increases in entropy occur in the blade row upstream 
of the stations where equations are solved and the entropy 
change along a streamline between blade rows is zero. 

S. The boundary layers on the turbine casing are not 
accounted for. 

The method of solution is as follows: 

1. Assume initial radial positions of the streamlines. 

2. Obtain the axial velocity distribution by solving 
the equation of motion at the stator outlet. The velocity 
distribution into the stator is assumed to be axial, and 
uniform 

3. Obtain stator loss coefficients. 

4. Check overall continuity and adjust the inlet Mach 
number aS necessary. 

5. Check the between-streamline continuity, and adjust 
streamline radial positions as necessary. 

6. Repeat this process for the rotor. 

7. Re-cycle all the above calculations, accounting for 


Streamline curvature, and repeat until convergence is reached. 
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C. METHOD OF SOLUTION 

The computer code written by Macchi was originally run 
on the IBM 360 computer. The program consisted of a deck of 
over 2000 program cards plus over 60 data cards. Since the 
deck could not be located, it was necessary to re-type the 
program from the listing in Macchi's paper. However, since 
the IBM 360 computer was soon to be replaced in the period 
in which the work was to be carried out, an alternate 
computer was sought. 

The HP-1000 series mini-computer located at TPL was 
selected for two reasons. First, the machine used FORTRAN 
as did Macchi's program. Secondly, it would be a benefit 
to TPL to have the program immediately available on the 
laboratory computer. 

The first steps were to analyze Macchi's program, in 
detail, and then to run it using his example input/output. 
In analyzing the program it became obvious that the computer 
program listing given in Ref. 2, was not the one used to 
obtain the listed output. Numerous discrepancies were found 
in the listing, some of which would have prevented the pro- 
gram from running; others would have caused incorrect results 
to be obtained. A listing of these discrepancies is con- 
tained in Appendix E. When the program was understood and 
flow-charted, it was keyed-in at the HP-1000 computer ter- 


minal. However, modifications were required to accomodate 
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the program within the mini-computer disc-based operating 


system. 


D. MODIFICATION TO THE COMPUTER CODE 

Since there was no card reader, variable input data such 
as turbine speed had to be entered using data or specification 
statements. This contributed in part to the most difficult 
problem, that of program size. The HP-1000 mini-computer uses 
a disc with a storage capability of 19.5 mega-bytes. However, 
the machine memory is only 124 K Bytes, of which only 29 K 
Bytes is available to a programmer. Also, the available 
memory is divided up, or partitioned into two 18 K and one 
11 K partitions, so that no single program can exceed 18 K. 
It was estimated that Macchi's program was over 100 K. So 
it was clear that the program would have to be modified if 
it were to run on the mini-computer. 

The first modification was to remove all subroutines 
from the program that were not actually used. It will be 
recalled that Macchi's program contained five methods for 
calculating loss coefficients and two methods for calculating 
gas outlet angles. It was decided that only the Traupel 
method of calculating loss coefficients would be retained. 
Traupel was selected for two reasons. Firstly, it was the 
method used by Macchi in his example calculations and there- 
fore the modifed program should still reproduce Macchi's 


results. Secondly, the method of Traupel is widely respected. 
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The method of calculating gas outlet angles was totally 
changed. Neither Ainley and Mathieson [Ref. 6] nor Traupel 
[Ref. 7] was used. Both methods required prohibitively 
large sections of computer code. The method selected was 
that of Vavra [Ref. 1]. 

Use of Vavra's method greatly simplified the program 
because this method predicts gas outlet angles independently 
of the inlet Mach number. Macchi's approach was to use 
Traupel's method which is dependent on the Mach number of 
the flow into the blade. 

The above simplifications reduced the program size from 
2287 lines to less than 1800 lines. However, this was still 
too large and the program could not be loaded without over- 
flowing the memory. 

The solution to the problem was found in program 
segmentation. In this process, the computer code is divided 
into a main program and several segments. Each segment is 
a "piece" of the original program. The segments are indi- 
vidually compiled and loaded. However, the segments are 
placed into memory only as they are needed to execute the 
overall program. Thus, a very large program can be made 
to run in the available 18 K partition. Since the present 
program was not originally intended for a mini-computer, 
segmentation was not straight forward. The method finally 


arrived at is detailed in Appendix C. Basically, the main 
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program consists of all the subroutines, while the three 
other segments contain coding which enables program flow to 
proceed in a logical manner. 

Successfully segmented, the program was run using Macchi's 
input. An output was obtained which agreed almost exactly 
with Macchi's results. All output quantities were within 1% 
of Macchi's quantities. The differences were, in all proba- 
bility, due to the different method of calculating gas outlet 
angles. 

After verifying Macchi's program, the drive turbine geo- 
metry was input and the program was run for a given set of 
operating conditions. The results are discussed in the follow- 
ing section. Note: The "verification" of Macchi's program 
amounted to verifying that the computer code now loaded into 
the HP-1000, was indeed Macchi's code. It was not known 
whether Macchi's output data were a good or bad prediction 


of performance since they were not compared with test results. 
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III. RESULTS OF AXIAL TURBINE PREDICTIONS 


A. USING BOTH COMPLETE AND MODIFIED PROGRAMS 

The drive turbine geometry was input and the following 
solution flow path was selected: 

1. Stator and rotor loss coefficients were functions of 
pressure ratio. 

2. The blockage factor, €*, used in the equation of con- 
tinuity was equal to the total loss coefficient. 

Four operating points were selected to test the validity 
of the program. Three were off-design points at which measured 
data were available and the fourth was the design point itself. 
Table II contains details of the selected test points for 
Run 1. 

The program variables were then changed and the following 
new solution flow path was selected: 

1. Stator and rotor loss coefficients were those calcu- 
lated by Traupel's method. 

2. The blockage factor, €*, was equal to the profile loss 
coefficient. 

After reviewing the results of Runs 1 and 2, a further 
modification was made to the program. The original program 
contained a subroutine which checked between-streamline con- 
tinuity. If the total mass flow rate at the stator and rotor 


exits was not evenly divided between the five streamlines, 
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the radial positions of the streamlines were adjusted and all 
steps were recalculated using the new streamline positions. 
Hence, for Run 3, a subroutine was removed and the main pro- 
gram was modified so that between-streamline continuity was 


not examined. 


B. COMPARISON WITH MEASURED DATA 

The results of Run 1, 2, and 3 are tabulated in Table III. 
Run 1 showed predictions of mass flowrate which departed about 
6% from the measured data. However, the horsepower predictions 
were off by as much as 16.17%. Furthermore, the computer pro- 
gram was unable to reach a solution for the design point. 

Run 2 produced worse results as is evident from the table. 
Again, the program was unable to converge to a solution at 
the design point. 

Run 3 produced more acceptable data. Additionally, con- 
vergence to a solution was noticeably faster and a solution 
was obtained at the design point. Because of this, the method 
used in Run 3 was used to map the drive turbine performance. 
The computer program used to obtain the results of Run 3 is 
described in detail in Appendix A and is listed in Appendix G. 
The results of Run 3 are shown plotted in Figures 3 through 8. 

To obtain the plots in Figures 3 and 6, a value of the 
total inlet temperature was approximated by the method of 
Vavra as contained in Ref. [14]. It was assumed that the 


static turbine discharge temperature should not be less than 
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45°F (505°R). This corresponds to the approximate temperature 
at which condensation of moisture in the air, assuming 100% re- 
lative humidity, will occur. The inlet temperature was given by 

Total Inlet Temperature = Static Outlet Temperature 

1- nytt -(h) Gy 
T i 

where Ns? the total-static turbine efficiency was assumed 
to be 81%, and Sp» Was the total to static pressure ratio. The 
total inlet temperature corresponding to each pressure ratio 
is given in Table IV. 

The computer output corresponding to each point on Figures 
3 through 8 is contained in Appendix F. Only one side of the 
dual flow turbine was analyzed, thus, the resulting printed 
values of horsepower, referred horsepower, moment, referred 
moment, flow rate and referred flow rate must be doubled to 
obtain the actual turbine characteristics which have been 


plotted in Figures 3 through 8. 
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IV. DISCUSSION 


The agreement of both the predicted flow rateand the horse- 
power vubtained in Run 3 with turbine test data was encouraging. 
It is to be noted however, that this agreement was obtained 
using a procedure which was conceptually incorrect. In Runs 
1 and 2, between-streamline continuity was checked and the 
streamlines were adjusted as necessary. In Run 3, between- 
streamline continuity was not checked, and as a result, the 
mass flow rate between streamlines was not precisely 25% of 
the total flowrate. It is noted however, that the deviations 
were less than 10.0% and while the radial positions of the 
streamlines varied by 10.%, the differences between predicted 
and measured output horsepower decreased from 24% to 4.5%. 
Since the enthalpy change on each streamline was computed 
using Euler's turbine equation, the total horsepower obtained 
by integration is sensitive to the streamline radial positions. 
On the other hand, the calculation of the overall mass flow: 
rate is primarily a function of the blade throat openings and 
inlet conditions of the flow. Consequently, in relaxing the 
requirement for between-streamline continuity, the output 
horsepower was changed significantly, while the overall 


flow rate was not. 
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Using this procedure, which preserves overall continuity, 
a performance map for the turbine was produced (Fig. 3-8) 
which agreed well with the off-design performance measurements 
made at lower speeds (Table III). It is noted however, that 
the inability of the program in its original form to predict 
the measured turbine performance is not explained, and both 
the program itself and the data input for the geometry should 
be closely re-examined. 

The difficulty in obtaining convergence to a solution at 
some operating points above the pressure ratio of 2.0 is 
likely to be the result of choking occurring on one or more 
of the streamlines. This was suspected but not fully explored. 

Finally, although the program was eventually made to run on 
the mini-computer, the time required to put the program into its 
final form was excessive since the original program was not 
written with segmentation in mind. When the segmented pro- 
gram was completed, only one operating point per run could be 
obtained. Thus, excessive time was spent compiling and loading 
the program. The execution time for the program averaged 2 
minutes at the lower pressure ratios and up to 30 minutes at 


the higher ones. This would be unacceptable if many points 


were to be examined. 
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V. CONCLUSIONS 


The program for calculating the performance of a single 
stage axial turbine reported by Macchi was revised, corrected 
and segmented and made to run on the Laboratory mini-computer. 
When applied to the geometry of the air-drive turbine of the 
compressor test rig, selecting specific options for the repre- 
sentation of loss coefficients, the revised program failed to 
converge when design-point test conditions where input. Also, 
the computed horsepower was in error by as much as 24% when 
the program predictions were compared with specific test data 
obtained from the rig at off-design (lower speed) conditions. 
The revised program did however closely reproduce the results 
given by Macchi in his original report for a specific turbine 
geometry. 

When the requirement that the computed stream surfaces 
be such that they divided the flow exactly into equal 25% 
increments was removed, the program converged satisfactorily 
for design point conditions and gave agreement with test data 
to within 5% in flow rate and horesepower at off-design 
conditions. 

The complete performance map for the air drive turbine 
was obtained with the program following this revision. Based 
on the favorable comparison with data so far obtained, the 
map is likely to describe the per: »rmance to better than a 10% 


uncertainty. This is considered to be satisfactory for 
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sizing the turbocharger for the proposed compressor riz 
modification. 

The following recommendations are made concerning further 

application or development of the computer program: 

1. The failure of the program to converge before the 
final revision was made should be analysed closely, 
and the final revision removed if possible. 

2. The geometrical input for the air drive turbine (which 
was taken from drawings) should be reexamined and the 
physical dimensions of the blade rows themselves should 
be measured. 

3. Consideration should be given to putting the corrected 
original version of the program onto the IBM 370 com- 
puter so that, when successfully operating, a turbine 


map can be calculated with a single load. 
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TABL I 


TURBINE GEOMETRY 
(see Figure 2; Dimensions in inches) 


STATOR: 

Hub Radius 2.764 
Mean Radius 3.196 
Tip Radius 3.627 — 
Blade Chord 1.003 
Blade Suction Side Radius of Curvature 2.8065 
Maximum Blade Thickness 2252 
T.E. Projected Thickness .03 
T.E. Normal Thickness .0186 
ROTOR: 

Hub Radius 2.693 
Mean Radius 55205 
Tip Radius 3.837 
Blade Chord 1.003 
Blade Suction Side Radius of Curvature 2.8065 
Maximum Blade Thickness 2252 
T.E. Projected Thickness .03 
T.E. Normal Thickness .0186 
Tip Clearance .O1l(estimated) 


TABLE II 


MEASURED/DESIGN DATA USED TO VERIFY THE PROGRAM 


POINT _RPM “Tn: (ny our (Ry TO. (PST) ER: M(gen) HAP. 
1 18310 579.2 550.8 23.56 1.602 5.542 110.1 
2 15200 557.4 517.8 20.43 1.390 4.698 63.27 
3 21300 578.9 506.8 27.13 1.846 7.033 172.0 
4e 30500 640.0 sue VavaG. 28 10.85 485 


*Design Point 
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TABLE III 


COMPARISON OF PREDICTED TURBINE PERFORMANCE 


VS MEASURED PERFORMANCE 


RUN I 
HORSEPOWER 
MEAS. $%DIFF. PREDICT. MEAS SDIFF. 
5.542 6.09 99.5 110.1 9.63 
4.698 0.89 52.5 63.27 16.17 
7.033 0.009 163.64 172.0 4.86 
10.85 “7° N.C 485 se 
RUN 2 
5.542 9.35 90.92 11052 17.4 
4.698 4.29 49.76 63.627 -21.355 
7.033 3.80 130.76 172.0 23.97 
10.85 == N.C. 485 a 
RUN 3 
5.542 5.01 LES .Az 110.1 2.74 
4.698 0.81 61.96 63.27 2.09 
7.0335 0.10 179.68 172.0 4.47 
10.85 4.15 444.18 485 8.42 


NC: Computer program would not converge 
a large number of iterations. 


to a solution after 


TABLE IV 


VALUES OF ASSUMED TOTAL INLET TEMPERATURE FOR EACH 
PRESSURE RATIO. GIVEN IN FIGS. 3, 5, 6, AND 7 


— 


PRESSURE RATIO TOTAL INLET TEMPERATURE (°R) 
1.4 545.5 
1.6 562.6 
1.8 $77.3 
2.0 591.0 
252 603.6 
2.4 615.3 
2.6 626.1 
2.8 636.6 
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APPENDIX: A 


DESCRIPTION OF THE COMPUTER PROGRAM 


A-1. INTRODUCTION 

To enable the program to run on the laboratory computer, 
the program was divided into 4 parts; a main program and 3 
segments. A detailed discussion of program segmentation on 
the HP-1000 computer series is contained in Appendix C. In 
the description which follows, the program is treated as if 
it were one large program with many subroutines. 

The description follows the individual steps from 
start to finish in the analysis. A program flowchart is 
given in Figure A-1 and the FORTRAN symbols used in the 
program are listed in Tables A-I to A-IX. 
A-2. DESCRIPTION 

A-2.1 Input Data 

There are 4 basic categories of input data; tur- 

bine geometry, operating conditions, special data and pro- 
gram control parameters. Since there was no card reader 
input device on the computer, all datawere entered using 
either data or specification statements. Explanations of 
the turbine geometry, operating conditions, special data 
and program control parameters are found in Table A-I through 
A-V. The nomenclature for the blading is given in Figure 


A-2. 
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A-2.2 Initial Geometric Calculations 

The first calculation performed is to establish 
the S streamline locations at the stator inlet (station 0). 
The streamlines are initially positioned such that there are 
equal areas (25% of the total flow area) between them. Next, 
blade heights of the stator and rotor are calculated using 
the hub and tip radii of each blade. Blade spacings for the 
stator and rotor are computed at 3 Streamlines; hub, mean 
and tip. The blade spacing on the mean streamline for the 


stator is given by 


2 
S = z Rm (A-1) 
where S = Blade spacing 


Number of stator blades 


N 
ul 


Rm = Mean stator radius 
A-2.3 Calculation of Gas Outlet Angles 
Subroutine VAVRA calculates gas outlet angles 
for both stator and rotor. The method is that of M.H. Vavra 
[Ref. 1]. The equation programmed in the subroutine is 


~i f GQ 
xn coe [$4 (|-G))),, 


where ®@ = Gas outlet angle 
a = Throat opening 
S = Blade spacing 


t,= Projected trailing edge thickness 


This method is much simpler than that used by Macchi since 
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there is no variation in outlet angle with Mach number (for 
sub-sonic conditions). Therefore, once calculated, the stator 
and rotor exit angles remain unchanged. Subroutine VAVRA 
computes exit angles for the hub, mean and tip streamlines. 
The outlet angles at streamlines two and four are computed 
later in the subroutines STATR and ROTO2. 

Before printing the input data, the program 
calculates the mean throat opening for the stator and for the 
rotor. The ten equally spaced radii and corresponding throat 
openings (part of the input geometry) are fitted with a fourth 
order Chebyschev polynomial. A throat opening corresponding 
to the mean radius is thus obtained. In the present appli- 
cation of the program to the drive turbine, the mean throat 
opening was obtained from the design drawing of the blading 
shown in Figure 2. It was assumed that the throat opening 
varied linearly with radial position and hence the throat 
openings at other radii could be calculated. The resulting 
throat openings are shown in the computer output under the 
heading of "Input Prints". The design values of the stator 
and rotor throat areas were obtained from the original 
design notes of M.H. Vavra. 

A-2.4 Calculation of the FlowRate 

Subroutine CHAN is called to calculate the mass 

flow rate entering the stator. The equations used are as 


follows: 
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eo 


Tro 


Te eee : 
}+ He, voy 
Ve 4 Ge ORT (A-4) 
Pp 
P= = 7 (A-5) 
} + aa Me 
A-6 
P= P/RT cae) 
2 2 
Az [ Rie - Rives | (A-7) 
m= PAV (A-8) 


. To c 
M ef is the reference Aue wicnieess flowrate 


and is used to check overall continuity later in the program. 


A-2.5 Solution of the Equation of Motion for the Stator 


Subroutine STATR is called to solve the equation 
of motion for the stator outlet conditions. The equation of 


motion which is programmed is as follows: 


2 
AVEO) 205 cae he CR y*+ (Sr) 
AE t . -cos't (K 2rm SE) - Gereas 
. | 


Ci cos oH _ | GiHeosle 2 ay | doit 


- $I 
Yi? Vai ax: Yi Vow, AY, (A-10) 


where C, = Ae J (a constant to convert H, the enthalpy 


from BTU FT? 
tem 19 see) 


7 @\ = Axial velocity at_a streamline 
Vea (3) Axial velocity at mean streamline 


x = R C2) . Streamline radius 
! Rm Mean streamline radius 
2 Ven 
ds” ol ye | - CiH costes 
aX, i ay l- ¥i7 Va. 


Ci H costa (1- ) 
— = Stator loss coefficient 
(which is initially assigned an 
estimated value) 
The derivation of this form of the equation of motion is 
given in Appendix B. However, at this stage of the analysis, 
the streamline curvature is assumed to be zero. Therefore, 


the equation of motion becomes: 


d (Ln be -_ olols 2 a. C cost 
ps ds Pe = ~ATanedk, S ~ SS sinra, + ee eCe Ais 
d Xi AX, Xi Y “Ver 
e 1m 
dH : a C,H Cos od: ds* 
x 7 |! YN? Vor 1 dX a 


The equation of motion is solved when the value 
of ¥, at each streamline satisfies the equation. The solu- 
tion is to first put the equation in the form: 


d (In. *) = us Cx) 


AX, 


where I(X) consists of the right hand side of equation (A-11). 


(A-12) 


Integrating equation (A-12) yields; 
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a) 


4 


InN, = j I) dX, +n c* (A-13) 


Xo 


where 2£nc* is the constant of integration when x = 1 and i 


1. With these boundary conditions Eq. (A-13) gives 
! 


An cts - ( Li dX, (A-14) 
XO 
using Eq. (A-14) in Eq. (A-13), 


x i 
In Y. = ( EOda xs | Todd f (A-15) 


Xo 


which becomes 
% 


a 2 
Se ae: ( I(x) dX ae: 
t 
Taking the inverse natural log and the square root of both 


. xX 
sides ye ( Tindx: 
a /, 
S tiaee = (A-17) 


Equation (A-17) is the form of the equation of motion solved 
in subroutine STATR. Solution of the equation gives five 
values of Y, and thus the value of the axial velocity at 
each of the five streamlines. Initially, the value of Y 
is taken to be 1 and the value of S— is taken to be zero. 
In succeeding iterations, the calculated value of y is used 
to obtain a new value of oa and so on. 

After calculating five values of Yo the ‘stator 


exit conditions are calculated at each streamline from the 


geometry of the velocity diagram. The convention for positive 
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and negative angles and velocities is defined in Figure A-3. 


The required relations are the following: 


Va, = Va. (3)- Yi 


(A-18) 
Vin, = Var: Taner ee 
V, = Va J Os Ki a36) 
Ve. = -Va, [AR /2L | (A-21) 


where L is the axial distance between stations and AR is the 
change in radial position of the streamline. Vey the radial 
component of velocity, is taken to be zero at this stage in 


the calculation. 


a 
V; 4 Vv; + Ve 
V2 (A-22) 
7, 2 Veers 
’ 29 ICp (A-23) 
Tis 2 Tro - Tro-Ti 
Pa a es (A-24) 
PR.= P / fue ie (A-25) 
a Lay 3 ng 
P, Pro | a (A- 26) 


Miz Vi [fe geRT 27) 


After the above quantities have been calculated at each 


streamline, subroutine STATR returns to the main program. 
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A-2.6 Calculation of the Stator Loss Coefficients 
The calculation of the stator loss coefficients 
at each streamline is accomplished by subroutine ALOS1. 
The method of solution to obtain these loss co- 
efficients is that formulated by Traupel [Ref. 7}. In 
Traupel's method, the value of the total loss coefficient is 


given by 


a ee te (A-28) 
total profile wall remaining 


The calculation of € requires 9 subroutines. Figure A-4 


total 
describes the connection between the subroutines and subroutine 
ALOSI. 

The first step is to obtain the value of the 
total profile loss coefficient, ED" ED is defined by Traupel 


to be 


. > fo XmX 5 bs fa ‘i ., (A-29) 


where Eno = initial value of the profile loss coefficient 


a mach number correction factor 
x, = trailing edge thickness correction factor 
En = loss coefficient due to mixing losses and separa- 


tion losses 
Ee = loss coefficient due to fan losses 
The total profile loss coefficient is calculated 
in the following manner. First, data for initial profile loss 
(En) as a function of gas outlet angle (a,) for various values 


of gas inlet angle (a,) is read from an array (Fig. A-5). 
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This is done by subroutine TRAU] and functions XPO and YC. 
The values of Epo are contained in two arrays XPO1 (5, 8) and 
XPO2 (6, 8). This is because the data shown plotted in Fig. 
A-5 has been divided into two sets. One set is for values 

of a, between 40° and 80°. The other is for values of a, 
between 80° and 170°. The FORTRAN symbols for the two ranges 
of values of a; are ALFOI1(I) and ALFO2(I) respectively. The 
FORTRAN symbol for the gas inlet angle is ALFl (J) once the 
data pints selected from the plot are entered, fifth and sixth 
degree polynomials respectively are fitted through the data 
points. The value of Epo can then be determined for given 
values of a; and ag. 

The mach number correction, Xn is obtained from 
Fig. f£-5. Subroutine CSIM calculates the value of Xn using 
straight line approximations of the plot. 

Subroutine CID calculates the remaining terms in 
the expression for Ep: These are Xo» em’ ee: They are ob- 
tained from the data in Fig. A-6 using the linear interpolation. 
The abscissa of the curves for X, and em is either f or 1-f 


where f is defined as S 


fei- —*—— 
tC SING, (A-30) 
where 6 = normal trailing edge thickness. 
t = blade spacing. 
a1 * gas outlet angle. 
The loss coefficient due to wall friction, ew? 


is calculated using | 
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a 2 ee Sape 


where t - blade maximum thickness 


(A-31) 


2 = blade height 


This equation is programmed in subroutine CSIW. 


The value of ER is obtained using subroutine CSIR. 


ER is defined by Traupel to be an all-inclusive loss coeffi- 
cient which accounts for any remaining losses not previously 


defined. It is written as 


: =X. fo (A-32) 


ERO is an initial value of ER which depends on the value of 
@, where » is given by 
Vi SINA, 
> =. U (A- 33) 
in which v, = true velocity of gas 
v = blade speed 
A plot of ERO vs ¢ is shown in Fig. A-7. The correction X 
is a function of s/% where 
s = chord length 
2 = blade height 
and is obtained using the data in the lower half of Fig. A-7. 
The total stator loss coefficient is computed 
for 3 streamlines; those at the hub, mean and tip. 
The loss coefficients at streamlines 2 and 4 are 


obtained by linear interpolation. 
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A refinement to the stator loss coefficient may 
be applied depending on the input value of one program control 


parameter. The following 3 variations of bs are available: 


ta 
i+ f, J 
£ | ~ | 
¢ © Pro 
S. 7 | at 
re 
oz | 
Pro (A-34) 
: 
e > é, (A-35) 
and xt 
vw 
| + ts | 
1+ Jo 6* 
. = ) ze 
rr ee | (A-36) 


where € = loss coefficient calculated using the method of 
9 


Traupel y-! 
a 


a 
¥ y-! 
BY = i+ +— (.9) 
(A-37) 
The values of the program control paramenter 


required to select between options are given in Table A-V. 


— 


Before returning to the main program, subroutine 
ALOS1 calculates a value of &* which is a blockage factor to 
be used in the equation of continuity. There are three ways 


to define §*; they are as follows: 


Ce 
a o (A- 38) 


f : 3 
~ ie (A-39) 
% 


f 3 t (A-40) 


A-2.7 Solution of the Continuity Equation After Re- 
turning to the Main Program 


The overall continuity at the stator exit is 
checked. Subroutine FLOWR performs this task. The flow chart 
for FLOWR is given in Fig. A-8. In FLOWR the mass flow rate 
required by continuity is checked against the calculated mass 
flow rate. If the calculated flow rate does not agree with 
that required by continuity, adjustments are made to the axial 
velocity and/or the inlet Mach number, as will be explained. 


The mass flow required by continuity is 


Mm oe M reer 
REQD ~ 
Zs Am Rim 


(A-41) 


where MpEF = reference mass flowrate as computed in subroutine 


CHAN 


53 


2. = # of stator blades 
AS = mean stator throat opening 
Rn = mean stator radius 


The mass flow rate at each streamline com- 


puted in this subroutine is 


: Pre Tre ACD 
. |_'tre. Ate z $ 


Ac 
where Z is an area reduction coefficient defined by 
H KK 
Le +x Py 
Ho -1 + § (A-43) 


Z gives the percentage of flow area hketween the blades over- 
which it is permissable to assume a uniform velocity. The 
boundary layer on both sides of the flow limits the available 
flow area ani the backage facotr, Z. accounts for this. 
Equation A-43, Z is seen to be a function of the energy para- 
meter H*** and €*. €&* is the value of the loss coefficient 
returned from subroutine ALOS1. The energy parameter is 


defined as 
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63 Energy thickness 
Ree = =~ = = 
H Displacement thickness (A548 4) 


H*** can be written as 


; ¥ ves 3 4 
+ Eas She, Ney NE 
piace net 3mti Smal Imei  Tmel mt) 
ee a a A eee a ee ee 
2 3 4 
1 + _! + Xe Bf Xe 4 Xe ae Xe 
Xe-! met = 3mtl = Smt Tmt! Im 
(A-45) 
where: 
m= .15 
X = 1-Piy for unchoked flow 
E Pro 
h deat 
X= 1 (Popyy] ¥ for choked flow 
and : Y_ 
nll ee ee 
Perit ~ ! ye! 


The derivation of Z and H*** os given in Appendix B. 


The expression for ¢, the flow function, 


for unchoked flow is 


(A-46) 


(A-47) 
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After calculating for each streamline, the flow 
rate is integrated from hub to tip and the resulting value is 
compared with Meeqd’ If the two values of flow rate agree to 
within a specified tolerance (see Table A-IV) continuity is 
considered to be satisfied. Then, after calculating the total 
percentage of mass flow between adjacent streamlines, sub- 
routine FLOWR returns to the main program. 

If the flow rates are not within tolerance the 
program checks to see if the actual mass flow is to high. If 
it is to high, the value of the axial velocity is lowered 
proportionally to the difference between the actual and 
required flow rates. 

If the actual flow rate is too low, the procedure 
is more complicated. First, the flow is checked to determine 
whether choking has occurred. Streamlines one and five are 
checked. If the flow is in fact choked at those streamlines, 
the inlet Mach number is lowered and the program loops back 
to recompute the reference mass flow rate and repeat the 
complete procedure. 

If the flow is not choked, the axial velocity is 
raised proportionally to the difference between actual and 
required flow rates and subroutine FLOWR returns to the main 
program. 

A-2.8 Calculation of the Rotor Inlet Conditions 

Continuity having been satisfied through the 
stator, the rotor relative inlet conditions are calculated. 


In subroutine ROTO1, the following expressions are used: 
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Us = Bed R STATOR. 
1a R RoTOR 
War, ad Vw, rh U 
-} 
6: = Taw [as] 
W, Pa Va, 


MF 4 Ve, + Ww," 
T. sc ( T, + W, ) 
TE 2 ge Ter. Use U, 
2AJeTep 
= a 
Pais ce Tre | -l! 
TE - ! 7; 
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(A- 


(A= 


-$1} 


-$3) 


55) 


He = ( Tre )(-24) 


Where Trp? Por and He are equivalent temperature, pressure 


(A-56) 


and enthalpy respectively. 
A-2.9 Calculation of the Rotor Exit Conditions 
Calculation of the rotor exit properties follows 
the same procedure as was used to compute the stator exit 
properties. The process is outlined here with notable 
differences explained. Subroutine ROTO2 calculates the 
rotor exit properties. A flowchart of ROTOZ is given in 
Fig. A-9. 
The first step in ROTO2 is to solve the equat:on 
of motion for each streamline. The equation of motion in 


terms of relative quantities is 


d (Ln Ya") _ 


As‘ dga_ ss 2 2 4Um CoS 6 


2 Un Ua cos” Ba C, cos” Ba dHe 


2 a + 
Ys Vex Y,° Va, Xa 
| - _Si He cos’ Ba : _dsi (A-57) 


Ya" Vary d Xz 
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At this point in the calculation, streamline curvature is 


neglected. Hence, Eq. (A-57) reduces to 


d(LnYa) —_ dé. 2 P 
coat oi a = "dt Tan Ba aN xT NBs 


4 Um cos Ba Siw $y _ _2Um Uz cos" Ba 


Ya Ve, Ya” Vor, * 


£icos*6x_ olHe +| p- —CHe cos®ba | Se 
Y," Va, 4X2 y* Va, AX, 


(A-58) 


The derivation of Eq. (A-S7) is contained in Appendix B. 


Equation (A-55) is similar in form to Eq. (A-10). Hence, the 


method of solution is identical to that employed by subroutine 
' STATR. However, after solving the equation, the value of Y, 
| at each streamline is examined to determine whether or not it 
| falls into the range -2<Y¥ <2,0, Values of x greater than 

| 2.0 are set equal to 2.0 while those less than .2 are set 


equal to .2. Successive values of Y, at each streamline are 


compared, and when the values of successive iterations are 
within a specified tolerance (see Table A-IV), the iteration 


ends. The values of Y, are used to calculate the rotor exit 


conditions using the following equations: 


Vag = Vaa (3) %G 


| 

| (A-59) 
| 59 

\ 

| 


Cos Ba (A-60) 
We (- Var ) D-CL 
« 5. (A-61) 
2 
Ta = Te -_ Wa 
Ae J cp (A-62) 


ae 
g 
pp 
\ 
= 
ry 
>| 
z 
“T 
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(A-63) 
Was= Vu, + U 
oe (A-64) 
Tas = Tre - Tre-Ta 
= (A-65) 
I- Fe 
a 
t-l 
Pee Pe ls 
a TE Tre (A-66) 


Subroutine ROTO2, then returns to the main pro- 
gram. 
After calculating the rotor outlet conditions, 


the rotor loss coefficients are computed. Subroutine ALOS2 


en eC 
a ee ee 


calculates the rotor loss coefficients following the process 
used in subroutine ALOS1 for the stator losses. The principle 
exception is that a tip clearance loss is also calculated and 
added to the total loss coefficient. The tip clearance loss 
coefficient is obtained from subroutine ALEAK which uses a 
straight line approximation to the curve shown in Fig. A-10. 
Subroutine ALOS2 also computes values of &€* and one of the 
three refinements to ER: 

Subroutine FLOWR is called to check continuity 
at the rotor exit. If continuity is satisfied, the program 
continues. If not, the same procedure is followed as pre- 
viously explained for the stator outlet (Fig. A-1l). 

A-2.10 Accounting for Streamline Curvature 

All calculations to this point have neglected 
streamline curvature and assumed that the streamlines remain 
fixed through the stator and rotor (Fig. A-11). The radial 
shift in a streamline between stator inlet and rotor outlet 


can be written as 


AK= Geren —R RoTOR (A-67) 
INLET OUTLET 


This is the net radial shift in a streamline between stations 
‘O' (stator inlet) and '2' (rotor outlet). It is shown in 
Section 16.4 of Ref. [5] that the radial shift in a streamline 


between the stator and the rotor (station 1) can be written 


as 
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| me 
< Reh eae zy R stator R ROR | (A-68) 
OuTLET INLET OUTLET 


The angle between the meridional velocity Ve and the axial 


velocity vy is X}. The radial velocity Ne, can be expressed as 


Ve = Va Taw d 


(A-69) 
and from Fig. 16(1) of Ref.[5], it follows that 
\ - AR | 
Tan re aL (A-70) 
Using Eq. (A-68) in Eq. (A-67), 
y, AR 
Ve = ae a (A-71) 
where a = Average streamline slope 
Also, from using Eq. (A-68) 
gL 
Cos N = . 5 (A-72) 
AR® +(QL) 
Rearranging; . 
| 
(A-73) 


The remaining term used in the calculation of streamline cur- 


vature (Section 16-4 of Ref. [5]) is 


6R 
Kr 


where K is the so called curvature factor. It usually has 
a value between 4 and 6 and in the program its value is taken 
to be S. Having calculated cos\, AR and 6R, the program 
repeats the solution process. However, the only quantity 
which is unchanged is the reference mass flow rate Met: In 
subroutine STATR the equation of motion is solved, this time 
accounting for streamline curvature. The same is true in 
subroutine ROTO2. 

The flow path of the program is identical to the 
section which did not account for streamline curvature. Next, 
the program computes an average pressure ratio at the rotor 


outlet using the expression 


Po, = Ps. r oa P, 
Fro Fro / STREAMLINE 4 Pro} cu. 
2 
Bb) ,{f\ , /& 
Pro Pro Pro 
SL. S.L. S.L- 
(A-74) 


If this pressure ratio is within a specified tolerance to 


the actual pressure ratio (which is input data) the program 
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ear ao as 


proceeds to the final stage of the calculations. If not, the 
inlet mach number is adjusted by an amount which depends on 
the difference between the calculated and specified beescure 
ratios. If the calculated pressure ratio is too high, the 
Mach number is lowered using 


M =m. - Pressure Ratio Difference (A-75) 


If the computed pressure ratio is too low, the Mach number 
is raised using 


Mo=M + Pressure Ratio Difference 


" - (A-76) 


In both cases, the program loops back to subroutine CHAN and 
proceeds to compute a new reference mass flow rate based on 
the new value of the inlet Mach number. The entire process 
is then repeated until the pressure ratios agree within the 
specified tolerance. 
A-2.11 Final Calculations 
Stator and rotor outlet conditions not previously 


calculated are computed as follows: 


Xa 


u 
>! 
2 


Va (A-77) 


Va 


(A-78) 


V3 ras vA + We 
Ak a U Vian -Ua Vika 
qev 
Tes = Tro cm ah 
=P 


yp 
u" 
= 
pp 
eran | 
| 
oli 
Ly 
7] 
el 


oT 
Pr, — P, Tro Le 
T; 
TL 
Ais = Tro Fa e 
Pro 


ROTOR EXIT 


RELATIVE MACH # 7 
Ts y — 


(A- 


(A- 


-81) 


-84) 


-85) 


Tro - Tris 
N = Ito- Tt, 
T-S i See 
To-7 Ta iS 
Stator Blade Tere= Ti 
Efficienc : 
y Tro- Ths 
Seton Bees Toe sat 
Efficiency 
Tro -Tais 
»* ~— 
vr = ling = Ler 
Tro -Tais 


Head Coefficient = 24 T (Tre Tans) 


U 


-1 
Blade-Jet Ratio = [Head Coefficient] 


Stator Exit . W, 
et a es 


Relative Mach # 
¥ Rge Th 
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(A- 86) 


(A-87) 


(A- 88) 


(A-89) 


(A-90) 


(A-91) 


(A-92) 


(A-93) 


(A-94) 


The turbine horsepower is obtained by integration. 
The Ah term at each streamline is weighted by the percentage 
of mass flow at that streamline. The product is then inte- 
grated from hub to tip and result, Ah, is used in the tur- 


bine horsepower equation 
Ah: J: m 
en 


550 (A-95) 


HP. = 
The moment is calculated using 


(H.P. )( 550) 


OD (A-96) 


M = 


Referred horsepower, moment, mass flowand RPM are calculated 


using 


Her, 


‘REF o§ (A-97) 
Meer = os (A-98) 
m ReF = . 2 (A-99) 
RPM pce = a (A-100) 
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where "Tee 


G 
7) 
A, 
% 
a 


14,7 


The values of the total-static efficiency, total- 

total efficiency, total-static pressure ratio, total-total 

pressure ratio, head coefficient, blade/jet ratio, r* and 

inlet mach number are then averaged. : 
With all calculations completed, the results are 

printed under the heading "STATOR SOLUTION", '"tROTOR SOLUTION", 

and "OVERALL TURBINE CHARACTERISTICS". 
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TABLE _ A-I 


TURBINE GEOMETRIC INPUT DATA (STATOR) 
(see Figure A-2; Dimensions in inches) 


FORTRAN SYMBOL DESCRIPTION 

zs Number of blades 

RS (1) Hub radius at stator outlet 

RS (3) Mean radius at stator outlet 
RS (5) Tip radius at stator outlet 

C Blade chord (mean) 

CI Blade chord (hub) 

co Blade chord (tip) 

E Blade curvature (mean) 

EI Blade curvature (hub) 

EO Blade curvature (tip) 

T Maximum blade thickness (mean) 
TI Maximum blade thickness (hub) 
TO Maximum blade thickness (tip) 
TE Projected T.E. thickness (mean) 
TEI Projected T.E. thickness (hub) 
TEO Projected T.E. thickness (tip) 
™N Normal T.E. thickness (mean) 
TNI Normal T.E. thickness (hub) 
TNO Normal T.E. thickness (tip) 
Al1(1-10) Ten values of throat diameter at 


10 equally spaced radii 
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FORTRAN SYMBOL DESCRIPTION 

AL Blade camber line length (mean) 
ALI Blade camber line length (hub) 
ALO Blade camber line length (tip) 
RC(1) Hub radius at stator inlet 
RC(3) Mean radius at stator inlet 
RC(5) Tip radius at stator inlet 
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TABLE A-II 


TURBINE GEOMETRIC INPUT DATA (ROTOR) 
(see Figure A-2; Dimensions in inches) 


FORTRAN SYMBOL DESCRIPTION 

ZR Number of blades 

RR(1) Hub radius 

RR(3) Mean radius 

RR(5) Tip radius 

CR Blade chord (mean) 

CIR Blade chord (hub) 

COR Blade chord (tip) 

ER Blade curvature (mean) 

EIR Blade curvature (hub) 

EOR Blade curvature (tip) 

TR Maximum blade thickness (mean) 
TIR Maximum blade thickness (hub) 
TOR Maxiaum blade thickness (tip) 
TER Projected T.E. thickness (mean) 
TEIR Projected T.E. thickness (hub) 
TEOR Projected T.E. thickness (tip) 
TNR Normal T.E. thickness (mean) 
TNIR Normal T.E. thickness (hub) 
TNOR Normal T.E. thickness (tip) 
TIPC Tip clearance 
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FORTRAN SYMBOL DESCRIPTION 
A2(1-10) 10 values of throat diameter at 


10 equally spaced radii 


ALR Blade camber line length (mean) 
ALIR Blade camber line length (hub) 
ALOR Blade camber line length (tip) 
cV Axial distance between stations 
CK Curvature Factor 


Te 


TABLE A-IIT 


TURBINE OPERATING CONDITIONS (INPUT DATA) 


FORTRAN SYMBOL DESCRIPTION 


AMC Assumed inlet Mach number 


Assumed stator exeit Mach number 


AMS (absolute) 
Assumed stator exit Mach number 


AMR (relative 

PTO Total inlet pressure (Pro) 

TTO Total inlet temperature (‘ro) 

PR Total-static pressure ratio 

RPM Operating speed (RPM) 

VA1(3) Assumed axial velocity in stator 

VA2(3 Assumed axial velocity in rotor 
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FORTRAN SYMBOL 


TOL 1 


TOL 2 


TOL 3 


TOL 4 


TABLE A-IV 


SPECIAL INPUT DATA 


DESCRIPTION 


Toler mce for convergence of 
equation of continuity 


Tolerance for between-S.L. 
continuity (not used) 


Tolerance in pressure ratio 
convergence 


Tolerance in equation of motion 
convergence 


TABLE A-V 


PROGRAM CONTROL PARAMETERS 


FORTRAN SYMBOL POSSIBLE VALUE EFFECT/MEANING 
Prints results in sub- 
IND 1 routines STATR, FLOWR, 
ROTO2 
1 No printing in the above 
ICL 1 Rotor is shrouded 
1 Rotor not shrouded 
IcOZz 1 E = Fo 
6 — = € (Y Pressure Ratio) 
8 BO Enea 


~N 


gm EDROFILE 


Wr 


f= STOTAL 
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= mee 


TABLE A-VI 


FORTRAN SYMBOLS IN THE MAIN PROGRAM 


FORTRAN SYMBOLS DESCRIPTION 
BESP ge = (1 + rae -(.8) 2} 
OI Stator throat opening (hub) 
00 Stator throat opening (tip) 
OIR Rotor throat opening (hub) 
OOR Rotor throat opening (tir) 
fe) Stator throat opening (mean) 
OR Rotor throat epening (mean) 
ANGZI Stator gas outlet angle (hub) 
ANG20 Stator gas outlet angle (tip) 
BETAI Rotor gas outlet angle (hub) 
BETAZ Rotor gas outlet angle (tip) 


fs es eee eewemanee Sa oygay EEE 
G Grav. constant, 32.174 Ppe-seoe 


CJ 778.16 FT.LBF/BTU 
EXP1 Vy a1 
EXP2 Y ayy 
ERRE Gas constant, 53.3459 a 
EMME Molecular mass, 28.970 LBM/LB MOLE 
GAM vy, Ratio of specific heats 
ETAT Total-total efficiency 
ETAI Total-static efficiency 
ETAS Stator blade efficiency 
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FORTRAN SYMBOL DESCRIPTION 

ETAR Rotor blade efficiency 

RSTAR Theoretical degree of reaction 

ALOS Head coefficient 

BLAJE Blade/jet ratio 

DR1 Radial shift of steamlines 

AMW1 Stator exit relative Mach Number 

AMS1 Stator exit absolute Mach Number 

AMV 2 Rotor exit absolute Mach Number 

AMR2 Rotor exit relative Mach Number 

DELH A 

HP Horsepower 

AMOM Moment 

THETA i) 

DELTA § 

HP1 Referred H.P. 

AMOM1 Referred moment 

RPM1 Referred RPM 

WLBM1 Referred mass flow rate 

ETAS Average total-static efficiency 

BETA6 Average total-total ~ “ssure ratio 

ETA6 ~—~—~—sC Average total-t- 11 . “+ iciency 

AKIS5 Average head coefficient 

RSTARS ~=~—«CAverage theoretical degree of reaction 
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ee re ee 


( 
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TABLE A-VII 


FORTRAN SYMBOLS IN SUBROUTINE CHAN 


FORTRAN SYMBOLS DESCRIPTION 
TTO Tro? total temp. at station # 
AMC Inlet Mach number 
PTO 7 Pro» total pressure at station § 
RC (1) Streamline radii at station # 
WLBM M, required mass flow, pAV 
TC Static temperature 
vc Velocity 
PC Static pressure 
RHO oe, density of air 
WCHAN Meng? reference mass flow 
WPERO % of M at each streamline 
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7 een 


TABLE A-VIII 


FORTRAN SYMBOLS IN SUBROUTINE STATR 


FORTRAN SYMBOL DE SCRI PTION 

ALFA1 Stator gas outlet angle 

X Ratio of streamline radius/ méan 
radius 

AMS Mach Number at station 1 

T Static temperature 

P Static pressure 

vl Absolute velocity 

VA1 Axial velocity 

Y Ratio of axial velocity to mean 


axial velocity 
S Entropy 


DSDX Entropy gradient between streamlines 
vul Tangential velocity 

PRAT (Total-static pressure ratio) 

T1IS Tits 


DALF Ix 


RSF Mean stator radius 

cea ge ae oe Na ae > NE ee he : 
DELR Stator in - ‘rotor out 

ZETAPS § 


ZETAS § 
Ss 


VR1 Radial velocity 
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TABLE A-IX 


FORTRAN SYMBOLS IN SUBROUTINE TRAUZ 


FORTRAN SYMBOL DESCRIPTION 
CSIP Xp» correction to 
R Sha initial profile loss coefficient 
ANG1 Gas outlet angle 
ANG2 Gas inlet angle 
Rl Xw Mach No, correction 

R3 Sp» remaining loss coefficient 

f rn oi eee 
RPRO Sp» total profile loss coefficient 

CL Rotor tip clearance 

YCL Tip clearance loss coefficient 

| RTOT Total loss coefficient 

| T Blade spacing 

DEZ Normal trailing edge thickness 

HM Blade height 
sa i ay 

PSID | Ses ieee ceeeet ere due to fan 

PSIF Sm? loss coefficient due to mixing 


and separation. 


UM Tip speed 


rte ee 


TABLE A-X 


FORTRAN SYMBOLS IN SUBROUTINE FLOWR 


FORTRAN SYMBOL DESCRIPTION 

PRATGR Critical pressure ratio 

PHICR ®oRIT? critical flow function 

HSTAR H***, energy parameter 

XI Z, area reduction coefficient 

PHI ®, flow function (unchoked flow) 

ARAT Streamline throat DIA/mean throat DIA 
AS Mean stator throat diameter 

AR Mean rotor throat diameter 


M required to satisfy continuity 
M calculated 


% of M at each streamline 
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TABLE  A-XI 


FORTRAN SYMBOLS IN SUBROUTINE ROTO1 


FORTRAN SYMBOL DESCRIPTION 
OMEG. uw, wheel speed (RAD/sec) 
U wR, tator mean 
U2 es wroter mean 
stator mean 
WU1 Wu see figure A-3 
BETA] Bi, see figure A-3 
W1 Wi, see figure A-3 
TTE Equivalent temperature 
PTE Equivaient pressure 
HE Equivalent enthalpy 
ZETAR - §, rotor loss coefficient 
ZETAPR §_, rotor profile loss coefficient 
DHEDX Enthalpy gradient between streamlines 


DSDX Entropy gradient between streamlines 


TABLE  A-XII 


FORTRAN SYMBOLS IN SUBROUTINE ROTO2 


FORTRAN SYMBOL DESCRIPTION 

BETA2 82, see figure A-3 

DBETDX as between adjacent streamlines 

VA2 Va.” axial velocity 

W2 W2, see figure A-3 

CL Axial distance between stations 

WR2 Radial component of velocity 

WU2 Wuze , see figure A-3 

Vu2 Vuz2, see figure A-3 

AMR Relative Mach No,at rotor exit 

T2 Te 

T2S Toe 

P2 P, 

PRAT 2 [Total-static pressure ratio] 
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INPUT DATA AND APPROXIMATIONS 


SOLVE EQUATION OF MOTION AND COMPUTE 
STATOR EXIT PROPERTIES 


SOLVE EQUATION OF MOTION AND 
COMPUTE POTOR EXIT PROPERTIES 


CALCULATE STREAMLINE. 
CUPVATURF AND SLOFE 


RE-CYCLE ALL OF THE ABOVE 
ACCOUNTING POR STREAMLINE CURVATURE 


ADJUST M, beers : 


YES 
PRINT RESULTS 


END 


FIGURE A-1l: PROGRAM FLOWCHART 
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« CAMBER LINE LENGTH 
= CHORD 

2 THROAT DIAMETER 

# CURVATURE RADIUS 

= SPACING 

® MAXIMUM THICKNESS 


TEs TRAILING EDGE THICKNESS PROJECTED IN 
PERIPHERAL DIRECTION 


TN= TRAILING EDGE THICKNESS, NORMAL TO FLOW 
DIRECTION 


FIGURE A-2: BLADE NOMENCLATURE 
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FIGURE A-3: VELOCITY DIAGRAM NOMENCLATURE 
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io er Se 


MAIN 


PROGRAM 


FIGURE A-4: INTERCONNECTION OF THE SUBROUTINES IN THE 


TRAUPEL METHOD 
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FECHA RENT 
pole ANNAN 


Bese alanine a panes ne wea aia ty 


CORRECTION FROM TRAUPEL 
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| 

i 

! 

FIGURE A-5: INITIAL PROFILE LOSS COEFFICIENT AND MACH NUMBER 
| 

| 

‘ 


FIGURE A-6: 
| 
4 


a re RY 
U a ae a7 46 as 


Abb. 8.44. Korrekturfaktor 73 und Mischverlust (5, 
infolge endlicher Austrittskantendicke oder 
Abléeung 


T.E.THICKNESS CORRECTION FACTOR, MIXING LOSS 
COEFFICIENT AND FAN LOSS COEFFICIENT FROM 
TRAUPEL 
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= 


_————- ee ee EE — 


FIGURE A-7: 


Untere Werte nur bei Profilen, 
aie wenig empfindiich 
aut Zustromwinkelanderung 


Abb. 8.4.7 Randverlust in cabie ee 


"REMAINING" LOSS COEFFICIENT FROM TRAUPEL 
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sie 2 fas So, 2 SA 


CALCULATE CRITICAL PRESSURE 
RATIO AND FLOW FUNCTION F 


Do FOR ZACH S.L. 


{ 
xe | - ( Rese ee 


(0 = CRITICAL FLOW FUNCT-ON 


es. _——— - 


ADD TIP CLEARANCE 


FIGURE A-8: SUBROUTINE FLOWR FLOWCHART . 
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ne EEEEEEEAEREREEnEEEnneE 


CALCULATE ROTOR OUTLET ANGLES 
FOR STREAMLINES 2 AND 4 { 
CALCULATE 38 POR EACH STREAMLINE 
dX 
DO FOR EACH S.L. 


SOLVE THE EQUATION OF MOTION 
AT THE ROTOR EXIT (STATION 2) 


| = 
ee 
t 
| os ae 
SOLUTION 
CONV. 
, ? { 
: YES 
CALCULATE ROTOR EXIT CONDITIONS : 
RETURN 
{ 
FIGURE A-9: SUBROUTINE ROTO2 FLOWCHART 
92 ; 
le 
| 
rr 
| 
| Y 


Dip 


Che 


Abb. 8.4.11 Faktor Ky, fiir Spaltverlustberechnung 


FIGURE A-10: TIP LEAKAGE LOSS COEFFICIENT PLOT FROM TRAUPEL 
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APPENDIX: B 


DERIVATION OF EQUATIONS USED IN THE PROGRAM 
B-1. EQUATION OF MOTION FOR RELATIVE FLOW: 


The equation of motion for relative flow Ref. [5] is 
— —> —_ 
Ve» WX(VXWrdIW)+ TVS 
(B-1) 
Using cylindrical coordinates, the terms of EQN (B-1) may be 


expressed as follows: 


Ve = to. de, ¢, oHe 7) ote 


r 98 = a2 or (B-2) 
te C2 ty 
r r 
VxXW = cue Boo ce 
ge oz or 
r Waa We Wr 


is dWr 2 Wa, 1 te d(rWu)  OWr 
az or r or a8 


AD=A10} lu2 NAVAL POSTGRADUATE SCHOOL MONTEREY CA F/G 972 


COMPUTER EVALUATION OF THE ON-AND@OFF DESIGN PERFORMANCE OF aN \retctu 
MAR 81 R CIRONE 
UNCLASSIFIED a 


gee 
Souoe 


~~ 


A NE ee 
« 


wx (7 XW) = 


Wa 


oY 


ir o 4, (3 We) ae) _ We. (2M - ae) 


Ls. ee. 

Wa. Wr 
1 | 3GrWa) Wel | | We = (Wa) 
r or 0} v | 390 > 


Se — Beal) _y, 5 [MEN awe i 


a2 or BY 
aWa. (We  >(rw 2 } = 
es) Wan = (28 aon) 


ue 


(B-4) 


(B-5) 
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; ' Qs ; As . a5 
TVS = Tlie bass a +c] 
(B-6) 
Combining equations (B-{1) through (B-6) the terms in (B-2) 
can be written as: 


= d= 
(B-8) 
PHR Wa | 3(rWu) 2Wp 
or re or = YS oe 
aWr oaWe a§ 
W 
- | at or [+201 Waser — 
(B-9) 
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Since the flow has been assumed to be axisymmetric, all deri- 
vatives with respect to ® are zero. Thus, Equations (B-7), 
(B-8) and (B-9) reduce to, respectively: 


_ We dlrWe) Wr olrWe) _ 2 Wy 
O= r 2 r ar onal 


0 Gr Wau) + Bi ‘ OS (B-11) 


ay ae ge ea 
as 
wo Wa + Pia 
a Tv or (B-12) 
Equation (B-10) may be written as 
Cr We) Wr 9 CrWa) a wr Wr. 
ae = We or a (B-13) 
Substituting into equation (B-11), 
dhe We We alr Wa) yyy, Wr 
ee ro Wa or oe 
Wr 2% ayy WaeWe _ - 2S 
ary We v az (B-14) 
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Multiplying equation (B-9) by Wr and Equation (B-14) by Wa 


results in 


w, 2Ha = Wa We 2(rW) Wa Wy ar e 


or r or 
Wa Wr one 4 2W WaW, + Wr T a (B-15) 
and 

dHe —_ _WaWre alrWw) Q\Vr 
Wea = is = er aa + WaWr > 


Wo, Wr ee — 2WWeWre + Wa T 2 


Adding these two equations yields 


dHe our gs 
We Sot + Wa SE = T] Wr 3+ We | a 


Since the flow has been assumed to be adiabatic, the total 
relative enthalpy, Hp» is constant along a streamline. Thus, 


HR dHrR 
Vu SO 2 Wass tr Ve Ose 
He ie or (B-18) 
and re-arranging, 
aoHe - _ _Wr dHe 
a2 Wa or 
(B-19) 
From equation (B-19), eq. (B-17) can be written as 7 
as Wr as 
aos eee OU B-20 
az We ar 
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Substituting Eq.s (B-19) and (B-20) into equation (B-15) gives 


_ We dHe We Wr 2GrWe) | yy dWr _ 
Wa or ro We +r Co 32 
2Na Wa W W as 
W r : 
or a Wa Wa or 
(B-21) 


Wi, 
Multiplying Equation (B-21) by an a yields 
R 


od He —- We _o(rWu) _ yy oWe 


or r or aan Y - 
We aha t+ QWWeer T i (B-22) | 


This expression is identical to equation (B-21) and is the 


Se, a ee 


equation which must be solved. It must be put into a form 
which can be solved by the computer. Re-writing equation 


(B-22) given that 


aWe eee aCe) 
Wa BWa = La Cwe') 
| . 2 Wa”) Wr , 2Wae 2CrWw) 
yields =e ee - 2Wea “>> + ee = ae = 
oHe aS) ces 
WJ - = 
4 Wa ak dr + at sr O (3-23) 


The relative enthalpy can be written 


| 
w,* . 2 2 
Ha= het 2e Bet ager ~ “% ager Ter 
(B- 24) 
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re ee 


The equivalent enthalpy, defined in ref. [1] is 


z 2 = 
H : e \. * W, ah Uz oa U, 
1 24eT 290 T 


(B-25) 
Hence, the relative enthalpy can be written as 
ge 
— a 
He= He - > 
(B- 26) 


Also, the turbine outlet static temperature can be written 


as 
_ He Wa 
Ts ~ Ce ACp 


(B-27) 


Substituting Eq. (B-26) and Eq. (B-27) into Eq. (B-23) and 


applying Eq. (B-21) to the rotor exit, gives 


2 
9CWe") - 2Wa, ce +2 Warr o Un, Warr ) 4 


rae a Caos 
a 
Yow Way,-2 sane ~ Vr He — 
fe He 2 +2Q Cp 
Wa" 0 Sa =O 
ac ei 
J on (B-28) 
Given the relationships: 2 
2 Wr 
Tan” \= 7 


] 


LN ae 


and 


Equation (B-28) can be written as 


O(Was) ayy re — War, aS2 


aes 3 Cpe ha Oly 
Wag , 24% Nay) He 
ee Pars vee = = 
20a), | oy w?] 2S 
rae’ Cp 2. a ais) 
2 
and substituting oiehs aur, into equation (B-29) 
gives : 
(Was) —-2Wea aWri = Wan oS2 -_ 
| a ¥ az Cpcom)1 ate 


Wava _3Cr Wu) He _ 
2 ae ane eT eae Tai + Yw Waw3-2 sh + 


au" fat = E = hs | a V2 


(B-30) 


, T 
Multiplying Eq. (B-30) by (2) results in the dimension- 
W 


a 
m 


less form of Equation (B-29): 


Yam 3 (Wa) = 2 War Wor Y. PY (Wr) — 


Wax, on. Wer, War ™ 32 
: 2 Wate War 
Wor Yam Sz ie) Wawa, Waa. V2 2 Ltom Waray Wa. 
2 2 : 
Wasa Cp COS 2 on Wa.,,, Vae V2 3 Cre 7 r n) 
+4 haw ee Wes 25 Yom oHe - awh ham 
Wee am Var Worm oh. Waren 
Tam a He _ War a Was oS2 
—_— |} > ee 
Ce es. Wasen Wasr on 
(B-31) 


Introducing the non-dimensional quantities 


Y= Wa (B-32) 


(B-33) 


. S (B-34) 
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Equation (B-31) is written as 


acy”) <a Y’ Fee aWr = Y* i5* 
ox Wa a2 cos7 dh — sso BX 
ay TANS 9 (XY Tans) UmY¥Tan@  _ 
Hl aHe ra 5. oPe Uy aHe = 


2 awe |_os* 
Y rants | 3st = O 


(B-35) 


The fourth term of Eq. (B-35) is 


2Y mu, abort) . ay Tov | xy 


d Taw® oY 
> + X TaN 6 5X 


= i TANG _? Taw B + 2Y Taw 6 oY i 
oX ox 
a 
Y¥ 2 
a X Tan 6 


104 


also, 2 Tew B : ) 8 
ox s 
an 


" 2Y Tan? B ot = Tay 2@ acy) 


Therefore, equation (B-35) can be written 


20) (i+ Ton'g)- a Yr, Be _Y* ast 


Wa, io aes oye 
2 TewB 08 ., Yer a YUmY Tan B® _ 
aN eta ay te ye Wee mm 


2 aHe Um U2 2He 2 a,| as* 
HS Ee | SE _ yy? ine 
Wa, o¥ War, | WF Tig | OX 


(B-36) 
Multiplying Eq. (B-36) by ese ) and observing that 
(| + TAN? 8 = moss) 
2 
ae + cos | va owe a = a 
ATaANA 26 + =~ SIN°B + “Rm Si ces 8 + aUmUa cos B u cos. 
A Wa'im Y 
acos’® QHe aHe cosB ay | 35” 
seg aie ate |B © 
(B- 37) 
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To account for streamline curvature the following terms are 


introduced: 


ie 


(By 


cos? } = 


(B-38) 
where A, the angle between the axial and radial components of 
velocity at a point, is approximated as the average value be- 
tween two stations. 


Also, 


(B- 39) 
where 6r is the streamline shift throught the rotor defined 


ace 6r = TRoror — TROTOR 


QUTLET INLET (B- 40) 
Substituting Eqs. (B-38) and (B-39) into (B-37) yields 


2 2 a 
dLnY) _ -c623 -k tm 8) — (4 ger) a 


~ 2 Tan 6 ae = SING =, Som UE EOS an Sit cos 
Y 


x = 


2UmVacos7B | 2cos® de _| aHecos’s _ sna] 3S 
dx 


We, 7 WarnY* Ak | War? 
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To obtain a dimensionless equation, the term 
Gi = AGT 


is introduced into Eq. (B-41) giving 


A(InY") al-(x 2087) sr) in, L? +n) | ol s*| _ 
Ady | ee aX 
dB 2 ce] 4Um Sin cos 8 Urn, Ua cus" 
a TAN —~ —— SIN + 
C,cos'® dHe Ci He cos °8 ds* 


oe = = 2 
Wee TO dd oh ae 83 EE ek 


(B- 42) 
Equation (B-42) is the form of equation of motion used in the 


computer program. 


B-2. EQUATION OF MOTION FOR ABSOLUTE FLOW 


The equation of motion for absolute flow 


(B- 43) 


Differs from the equation of motion for relative flow 


VHe= Wx (7 xW +a )+TVS 
- (B-44) 
only by the term Wx2W which is the Coroiolis acceleration. 
To obtain the programmed form of the eauation of motion for 
the stator, the previous derivation is followed, but with 


U = 0, He becomes H, W becomes V, and 6 becomes a. 


B-3 THE AREA RESTRICTION FACTOR Z 
The condition at the outlet of a blade row with boundary 
layers on both sides of the flow channel is shown in Fig. B-l. 
The flow is considered to be turbulent within the boundary 
layer while, outside the layer, the velocity of the flow is 
the theoretical velocity. Assuming a power-law velocity pro- 


i file, the velocity may be written, 


ve » ae as 


| (B-45) 
| The mass flow rate exiting the blade row can be expressed as 
| ; ¢ z§ 

| Mm = Cru Vru CosXd | S- za - + 

| ; S COSA, Cos oJ 


| 2 ( Ane ol y 
, (B- 46) 
where Opp. and Ven represent the ideal conditions for an 


isentropic expansion through the blade row to the discharge 


: 
| 
| 
| 


pressure Pd, which is assumed to be constant across the blade 
spacing. The discharge angle of the flow leaving the blade 


row is closely approximated by the expression [Ref. 1] 


a OW. 
cose (B-47) 
Inserting Eq. (B- 47) into (B- 46) and reducing yields 
§ 2 
AL 
Me ze (ee 
™ Pru Vin a oO ( A (Pm Abny )) 
(B-48 ) 
Assuming a perfect gas 
P Ttx Tro — (Tre - Tr) 
ee LR) sae ae a eee 
pecs, tee = rae: 

Pr T To Co ak es) 
Defining y-| (B-49) 
(= r 

heels Pr (B-50) 
Equation (B- 49) can be written 
e |- Xe 
= = yn 2 
Fru |- Xe( ) (B-51 ) 
Vr 
Substituting Eq. (B- 51) into (B- 45) yields 
$ Cope 
Mh = Cun Of E(i-(-ney( = aa 
~ KE 


ce] 
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Using the displacement thickness given by 


g* 6. I~ (ive) (4 al h 


fo) (I-%e9"”) 


the mass flow ratecan be written as 


35* 


e . a _ 
mM Pra Vrw xy 


(B-54) 


The loss coefficient, expressed in terms of average kinetic 


energy lost is 


ee 


AE E 
f =m (vad) = 1 Vm 


a a 
(B-55) 
where E is the actual kinetic energy of the flow, given by 
V; sw 
fa 
E = PriV'ry(a- 25) Mt + 5 ( Pe > dy 
oO (B-56) 
Substituting Eq. (B-51) into (B-56) gives 
a | 3m 
E=@, da I- £1 U-xe)| —i___—d 
a. ) -Xen™”) 
(B-57) 
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The energy thickness is written as 


! 3m 
ek! ae 
CC te) ) (ore) 
(B-58) 


The loss coefficient can therefore be written as 
5 


f | fo Ze ee 
= S ¥ 
—_ ea Se B- 
l- 22 (B-59) 
The area restriction factor 2Z, is the fraction of the flow 
area through which the uniform theoretical velocity would 


produce the actual flow rate, thus 


>a ad 
~> CO 


B-60 
Defining the energy parameter (a form factor) as 
*¥e cs 
H = oe | 
(B-61) 
using Equations (B-59) and (B-61), Eq. (B-60) becomes 
py ¥ #4 7 | 
Zz? 
H wr ad +> 
(B-62) 


| 
where Ep is the profile loss coefficient. | 
| 
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B-4. THE ENERGY PARAMETER, H*** 
In Equations (B-53) and (B-58) the denominator of the 


integrand is expanded using the binomial theorem, so that 


an 7! etn 2 4 3 
(i-Xen )=1+Xen +YXe h o Xe noe 6-83) 


The integral of Equation (8-58) is now written as 


am d 3m om 2 7m 


oO 


Il im 


tet Xe A el AA (3-64) 


which, on integration becomes 


3m | X a 3 
(2 aly « piety Ne gike 


l-Xe ) 2m 3mtl Smt Tmt Amt! 
Ft XE | (B-65) 
Ilm+! ae 


swe |, Me ye Me 
Crna | - 3m+l Smt mH Gin +l 
4 
Xe (B- 66) 
Il mt (\-x,) 
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which can be written as 


a 
ge e -1) P l Xe Xe + 
& c Xe-l 3m+l Smt Tm +1 
3 4 
Xe , Xe 
Tm+l Il mel 
(B-67) 
In a similar Manner, 
Meg YE Xe? Xe° 
ae (%e-1) ———— + re se t 
§ -| M+ smt Sm+ Pmt 
ie wes 
Im +l 
(B-68) 
Substituting Eq. (B-67) and Eq. (B-68) into Eq. (B-61), the 
equation for H*** used in the computer program is obtained: 
a 3 4 
| Xe Xe RE Cg 
er Sma * Smet” Tmt” Amel mel 
ay 
a | a i a 
H 2 3 4 
| Xe Xe Xe Xe 


Xe- m+ 3mt Smt Tm 41 Am +1 


(B-69) 


FIGURE B-1: 


BOUNDARY LAYER EFFECTS AT THE EXIT 
OF A BLADE ROW 
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APPENDIX: C 
PROGRAM SEGMENTATION ON THE HP-1000 


Segmentation allows large programs to be run on the 
HP-1000. The program is divided by the programmer into a 
main program and several segments, which are stored on the 
disc. Each segment and the main program are then compiled 
and loaded. When the program is executed, the main program 
and its segments are called into memory individually, and 
only as they are needed for execution. In this manner, a 
program can run in a partition which is smaller than that 
program's total size. 

When the main program has performed all executable state- 
ments, the first segment is called into memory by an EXEC 
call. The system then loads that segment from the disc into 
a memory block following the end of the main program. The 
process is illustrated in Figure C-1. Note; the main program 
plus the largest segment may not together exceed 29 k. Once 

| a segment is in memory it can call another segment. 

| When executing, any segment can call any subroutine which 

: is attached to the main program. It was this feature which 

| allowed the present program to be run. All subroutines were 
consisted of nineteen subroutines and functions. A segment 


may not return to the main program. Communication of data 


115 


placed within the main program. In fact, the main program 
| 
| 


fay) ae ee 


between the main program and the segments is accomplished 
through a common block. 

The four segments of the present program are 'MAIN", 
"SHORT", "PART 2" and "PART 3". The manner in which control 
is passed from the main program to the first segment and from 
the first segment to the second is as follows: 


BLOCK DATA 


END 

PROGRAM THESS 

DIMENSION INAM (3) 

DATA INAM /2HSH, 2HOR, 2HT / 


CALL EXEC (8, INAM) 

END 

PROGRAM SHORT (5) 

DIMENSION INAM (3) 

DIMENSION NAME (3) 

DATA INAM /2HSH, 2HOR, 2HT / 
DATA NAME /2HPA, 2HRT, 2H2 / 


CALL EXEC (8, NAME) 
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END 
PROGRAM PART 2 (5) 

DIMENSION NAME (3) 

DIMENSION NAMR (3) 

DATA NAME /2HPA, 2HRT, 2H2 / 
DATA NAMR /2HPA, 2HRT, 2H3 / 


CALL EXEC (8, NAMR) 

END 

PROGRAM PART 3 (5) 

DIMENSION NAME (3) 

DIMENSION NAMR (3) 

DATA NAME /2HPA, 2HRT, 2H2 / 
DATA NAMR /2HPA, 2HRT, 2H3 / 


END 
The "(5S)" after the program name indicates that it is a 
program segment. Note the manner in which the program name 


is put into a data statement using the Hollerith notation. 
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SEGMENT 1 SEGMENT 
OVERLAY 


AREA 


saa 3 


MAIN 
PROGRAM 
AREA 


SEGMENT 
1 


SEGMENT 2 


PROGRAM 


PARTITION 
AREA 


DISC MEMORY 
SYSTEM 
TABLES 


LOGICAL MEMORY 


FIGURE C-1: PROGRAM SEGMENTATION-ILLUSTRATION OF THE MAIN 
PROGRAM CALLING A SEGMENT INTO LOGICAL MEMORY 
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APPENDIX: D 


RUNNING THE COMPUTER PROGRAM 


If the reader is unfamiliar with the HP-1000 Computer 
System, references [11] and [12] should be consulted before 


attempting to run the program. 


D-1. DATA INPUT 

Using the editor, input the following data into segment 
"SHORT". 

1. Turbine operating conditions: referring to Table 
A-III, type in appropriate data in lines 66 through 69 and 
74 through 78. 

2. Special input data/program control parameters: 
referring to Tables A-IV and A-V, type in appropriate data 
in lines 83 through 98. 

3. Turbine geometry: referring to Tables A-1 and 
A-II, type in data for stator and rotor in lines 103 through 


186. 


D-2 COMPILING THE PROGRAM 
1. To compile the main program type: 
:RU,FTN4,MAIN::25,-,- 


2. Compile the first segment: 


:>RU,FTN4,SHORT: :25,-,- 


3. Compile the second segment: 
:RU,FIN4, PART2::25,-,7 
4. Compile the final segment: 


RU, FIN4PART3: :25,-,- 


LOADING THE PROGRAM 
Type 
: RU, LOADR 
Tap return key 
Will display 
LOADR: 
Type 
OP,LB 
Will display 
LOADR: 
Type 
:RE,$MAIN::25 
Will display 
LOADR: 
Type 
:RE,$SHORT: :25 
Will display 
LOADR: 
Type 
:RE,$PART2::25 
Will display 
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LOADR: 


Type 
:RE,%PART3::25 | 


Will display 
LOADR: 
Type 
> END 


After the end statement, the loader will display that the 


program is ready for execution. 


D-4 RUNNING THE PROGRAM 


Type 
RUN, THESS 
The program will be executed and no further action by 
the operator is required. The computed pressure ratio of 
each iteration of the outer loop of the program is displayed 
on the screen as it is calculated. The operator therefore 
has some idea where in the iteration process the computer 


| 
program is executing. 
1 


— 


_ 
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APPENDIX: E 


DISCREPANCIES IN MACCHI'S PROGRAM 

1. Main program, lines 21 and 22; the value of ICL has 
not yet been read. 

2. Main program, lines 163-166; the Traupel method of 
calculating gas outlet angles does not take the Mach number 
into consideration. However, in lines 163-164, the program 
is attempting to draw a parabola through points which repre- 
sent outlet angle as a function of Mach number. 

3. Main program line 281; the calling of subroutine 
SLINE is questionable. Parameters are transferred to that 
subroutine, but many of them have not yet been defined 
(HE, DHEDX, WPER2Z, DSDX1). These undefined variables will 
be set equal to zero by the IBM 360 and 370 computers. Thus, 
in line 10 of subroutine SLINE, the value of DWDX will be 
zero and in line 17, division by zero will occur and the 
execution of the program should cease. 

4. Subroutine ROTORI lines 22 and 26; the stator radii 


are used in the calculation whereas the rotor radii should 


be used. 


5. Subroutine ASOSI, line 107; the correct Fortran code 
is 


ZETAPS(I) = .5 * ZETAS(I) 
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6. Subroutine ALOS2, line 121; the correct Fortran code 
is 
ZETAPR(I) = .5 * ZETAPR(I) 
7. Subroutine ALOS2, lines 123-126; the stator radii are 
used in the calculation whereas the rotor radii should be used. 
8. Subroutine ANGAIN, line 14; the correct Fortran code 
is 
AO = ATAN(1. -XCL/H*CH*COS (ANG1) /COS (ANG2) * 
TAN(ANG2) + XCL/H*CL*COS(ANG1)/COS (ANG2) * 
TAN (ANG1) 
Note: Since reference [2] was published, Professor Macchi's 
program has been further developed by Professor Macchi under 
private sponsorship [Ref. 13]. The new code however, is not 


generally available. 
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APPENDIX F 


COMPUTER OUTPUT 


ee 
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